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We  study  an  exactly  solvable  quantum  spin  model  of  Kitaev  type  on  the  kagome  lattice.  We  find  a  rich  phase 
diagram  which  includes  a  topological  (gapped)  chiral  spin  liquid  with  gapless  chiral  edge  states,  and  a  gapless 
chiral  spin  liquid  phase  with  a  spin  Fermi  surface.  The  ground  state  of  the  current  model  contains  an  odd  number 
of  electrons  per  unit  cell,  which  qualitatively  distinguishes  it  from  previously  studied  exactly  solvable  models 
with  a  spin  Fermi  surface.  Moreover,  we  show  that  the  spin  Fermi  surface  is  stable  against  weak  perturbations. 
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Introduction.  The  search  for  exotic  ground  states  of  quan¬ 
tum  many-body  systems  has  been  at  the  fore  in  experimental 
and  theoretical  research  for  a  number  of  decades,  yet  nature 
continues  to  tantalize  us  with  ever  more  interesting  behaviors. 
In  recent  years,  systems  with  some  type  of  topological  order 
beyond  Landau’s  broken  symmetry  principle,1-2  such  as  the 
fractional  quantum  Hall  effect3  and  quantum  spin  liquids,4 
have  garnered  a  great  deal  of  attention.  A  quantum  spin 
liquid  (QSL)  is  an  insulating  state  that  does  not  exhibit  any 
conventional  magnetic  order  at  zero  temperature,  typically  due 
to  a  delicate  balance  of  competing  and/or  “frustrating”  effects 
(for  a  review,  see  Refs.  5  and  6). 

Various  flavors  of  QSLs  with  different  types  of  quasiparticle 
excitations  and  braiding  statistics  have  been  proposed. 7~9  How¬ 
ever,  the  reliability  of  the  approximations  typically  invoked  in 
their  solution  are  often  imperfectly  understood.  Thus,  it  is 
useful  to  find  exactly  solvable  (even  though  often  contrived) 
models  with  spin  liquid  ground  states  so  as  to  establish 
their  stability  as  phases  of  matter  and  to  understand  their 
essential  physical  properties.  The  existence  of  such  states  is 
compellingly  established10  for  quantum  dimer  models1 1  on  the 
triangular  lattice. 

More  recently,  the  discovery  by  Kitaev12  of  an  exactly 
solvable,  interacting  two-dimensional  (2D)  spin- 1/2  model 
on  the  honeycomb  lattice  with  a  spin  liquid  phase  came  as  a 
great  boon  to  the  study  of  exotic  ground  states  in  magnetic 
systems.  Since  then,  variants13  of  Kitaev’s  model  on  trivalent 
lattices  have  appeared  which  realize  chiral  spin  liquids  with 
non- Abelian  anyons14  and  spin  liquids  with  a  spin  Fermi 
surface. 15-1 7  They  have  also  been  used  to  study  quantum 
critical  points,18,19  entanglement  entropy  and  entanglement 
spectrum,20  and  edge  solitons21  among  others.22-23  If  the  coor¬ 
dination  number  of  the  lattice  is  larger  than  3,  a  generalization 
of  the  Kitaev  model,  the  so-called  T -matrix  model.13-24  23  with 
extra  degrees  of  freedom  on  each  site  is  needed  to  ensure  exact 
solvability.  These  additional  degrees  of  freedom  can  be  either 
interpreted  as  an  “orbital”  degree  of  freedom,26'29  or  one  may 
consider  each  site  as  having  a  larger  spin,  but  no  orbital  degrees 
of  freedom.15,30 

Partly  motivated  by  the  recent  interest  in  finding  a  spin 
liquid  with  a  spinon  Fermi  surface  in  Herbertsmithite — a 
Mott  insulator  on  the  kagome  lattice31 — we  study  an  exactly 
solvable  spin-3/2  model  on  the  kagome  lattice  (see  Fig.  1). 


Our  main  finding  is  that  our  model  possesses  a  gapless 
spin  liquid  phase  with  a  finite  spin  Fermi  surface.  We  show 
that  the  spin  Fermi  surface  is  stable  against  any  weak 
perturbations.  The  stability  follows  from  the  lack  of  inversion 
symmetry,  n  rotation,  and  time -reversal  symmetry  (TRS)  (that 
is  spontaneously  broken)  in  the  ground  state.  Both  the  odd 
number  of  electrons  per  unit  cell  in  the  ground  state  and  the 
stability  of  the  spin  Fermi  surface  qualitatively  distinguishes 
our  model  from  other  known  exactly  solvable  models  with  a 
finite  spin  Fermi  surface. 15-1 7  We  also  find  a  gapped  Abelian 
chiral  spin  liquid  phase,14  which  possesses  two  gapless  chiral 
Majorana  edge  states  for  a  lattice  geometry  with  boundary. 

Model.  The  F-matrix  model  describing  S  —  3/2  spins15 
we  study  respects  the  translational  and  threefold  rotational 
symmetry  of  the  kagome  lattice  and  is  given  by 

h  =  jaJ2  r'r;  +  'v  £  r3r4  +  j5J2  r,5 

(ij)eA  (ij)eV  i 

+  J'a  £  r,15r-5  +  4  E  (D 

<0')e  A  (h')eV 


where  we  have  distinguished  the  nearest-neighbor  couplings 
J,j  as  JA,  Jy  and  JL  as  J'A  and  Jf  based  on  whether  the  link  (i  j) 
belongs  to  an  up  (A)  triangle  or  down  (V)  triangle,  and  (i j)  is 
taken  in  the  counterclockwise  sense  for  each  triangle.  Locally, 
the  five  T  matrices  satisfy  a  Clifford  algebra  {Ff  ,Ff }  =  2<5ab, 
where  a,  b  =  1.  ...  ,5,  and  rab  =  [Ta  ,Tb~\/(2i).  In  terms  of 
the  components  of  the  spin  5  =  3/2  operators,15-’0 


r1 

r4 


1 

71 

i 
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{5-v,52},  r2  = 

[(Yv)2  -  (S*)2], 


7!|S! 


r3  = 


l 

71 


{Sx,Sy}, 

(2) 


With  the  identification  (2),  it  is  clear  the  model  (1)  has  a 
global  Ising  spin  symmetry  under  180°  rotations  about  the  z 
axis,  and  possesses  TRS,  although  TRS  will  be  spontaneously 
broken  in  the  ground  state  as  we  describe  below,  in  addition 
to  the  translational  and  threefold  rotational  lattice  symmetry 
mentioned  above. 

Method  of  solution.  The  key  to  the  exact  solvability12  of 
the  model  (1)  is  an  infinite  number  of  conserved  operators  Wp 
that  satisfy  [WP,TC\  =  0  and  [Wp,  WP']  —  0  for  all  elementary 
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FIG.  1.  (a)  is  a  schematic  representation  of  the  kagome  lattice 
which  is  a  network  of  corner  sharing  triangles  in  two  dimensions 
with  three  sites  per  unit  cell,  (b)-(e)  are  four  energetically  distinct  flux 
configurations  preserving  the  translational  symmetry  of  the  lattice.  A 
doubled  (magnetic)  unit  cell  is  required  for  (d)  and  (e),  where  the 
total  flux  through  a  unit  cell  is  ±jr . 


plaquettes  labeled  by  p  and  /?'.  As  shown  in  Fig.  1,  the  kagome 
lattice  contains  elementary  up  and  down  triangular  plaquettes, 
and  hexagonal  ones.  Specifically,  Wa  =  r?2r]2r^2,  Wv  = 
r?4rf  r!4>  and  =  r23r]4rf  r/4r23r,;4,  where  the  sites 

are  taken  in  a  counterclockwise  fashion  and  for  Wq  the  first 
link  (i j)  is  assumed  to  lie  on  a  A. 

In  order  to  solve  the  Hamiltonian  (1)  we  introduce  a 
Majorana  representation  of  the  V  matrices,14 

r?  =  itfa,  r?  =  iCidi,  r f  =  igdif  (3) 

with  a  =  1,2, 3, 4.  There  are  thus  six  Majorana  species  on 
each  site  i :  {£I1,£2,fi3,§4,c;,<f1}.  The  Majorana  representation 
enlarges  the  spin-3/2  Hilbert  space,  so  that  one  must  enforce 
the  constraint  Z),  =  -r/r?r?rfr?  =  =  1. 

namely,  for  any  physical  state  i'I')phyS  AWphys  =  l^/phys 
for  any  i .  From  the  enlarged  Hilbert  space,  physical  states  are 
obtained  by  applying  the  projection  operator  P  = 
where  the  product  is  over  all  sites  in  the  lattice. 

Using  relations  (3),  the  Hamiltonian  (1)  becomes 

H  —  ^ [  Jp iuijCiCj  +  J'j iujjdjdj]  +  A  ^  ic,d,,  (4) 
<(/>  ■ 

where  u/j  =  — up ,  with  ipj  =  —iijfljj  if  ij  e  A  and  uij  = 
— if  ij  e  V.  The  original  Hamiltonian  is  obtained 
through  H  —  PHP.  Since  [«,,,?/]  =  [uij,Ui’y  =  0,  the  orig¬ 
inal  spin  Hamiltonian  has  been  reduced  to  a  model  of  free 
Majorana  fermions  moving  in  a  static  background  of  Z2  gauge 
fields: 1 2  In  (4)  the  utj  can  be  replaced  by  their  eigenvalues  ±  1 . 
It  remains  to  determine  the  pattern  of  the  Ujj  (up  to  gauge 
transformations)  that  yields  the  lowest-energy  state.  We  will 
limit  our  discussion  to  positive  J  couplings  only. 

In  terms  of  static  Z2  gauge  fields,  we  define  a  flux  <j>p  via 
exp (i(pp)  =  Y\  jkep  iujk,  where  jk  is  taken  counterclockwise, 
on  each  elementary  plaque tte  p.  It  is  clear  that  <pp  =  ±jt/2 
in  the  triangular  plaquettes,  and  <j)p  —  O.tt  in  the  hexag¬ 
onal  plaquettes.  Since  Wp  =  Y\(ij)ep  uij>  where  ij  is  also 
taken  counterclockwise,  exp (/</>p)  =  —  Wp  (— iWp )  when  p 
is  hexagon  (triangular)  plaquettes.  Under  TRS,  Wp  — >■  ±  W p , 
where  +(— )  is  for  hexagon  (triangle)  plaquettes;  it  follows 
that  tpp  — ►  —(f>p  for  triangle  plaquettes  while  tpp  remains 


unchanged  for  hexagon  plaquettes.  Consequently,  a  ground 
state  with  a  certain  flux  pattern  {<pp}  (regardless  of  any 
particular  choice  of  m^)  spontaneously  breaks  TRS  and  its 
energy  must  be  degenerate  with  the  ground  state  with  a  flux 
pattern  obtained  from  {<pp}  by  changing  tpp  —*■  —  (j>p  on  all 
triangular  plaquettes.14 

In  general,  the  nature  of  a  ground  state  and  its  excitations 
crucially  depends  on  the  flux  pattern  {<pp}  in  the  ground  state. 
As  we  will  see  below,  the  richness  of  this  model  is  due  in  part 
to  the  freedom  to  add  additional  terms  to  (1)  that  preserve  the 
exact  solvability,  but  can  be  used  to  select  different  (gauge- 
inequivalent)  configurations  of  w,-;-  as  the  ground  state. 

Chiral  spin  liquid.  A  chiral  spin  liquid  (CSL)  breaks  TRS 
spontaneously  and  has  no  conventional  magnetic  order.414  ’2 
As  we  have  discussed,  the  current  model  ( 1 )  indeed  breaks  TRS 
spontaneously.  It  remains  to  be  shown  that  a  ground  state  exists 
without  breaking  spin-rotational  or  translational  symmetry.  To 
determine  the  ground-state  flux  configuration,  we  compute  the 
ground-state  energy  Eg  of  the  Hamiltonian  (4)  as  a  function 
of  the  Ujj. 

For  a  bipartite  lattice,  the  ground-state  flux  pattern  is 
determined  by  the  Lieb  theorem”  and  the  flux  is  always 
uniform.  However,  the  kagome  lattice  is  not  bipartite,  so 
we  must  determine  it  numerically  and  analytically  in  some 
special  limits.14  By  diagonalizing  systems  of  up  to  ~104  sites 
and  computing  Ec,  we  indeed  find  that  for  various  possible 
parameters  in  Hamiltonian  (1),  the  ground  state  has  a  uniform 
flux  pattern.  There  are  four  possible  time-reversal-inequivalent 
uniform  flux  configuration  labeled  by  {(pA,<Pv  ,<Po)  = 
{f,f,7r},{  — f,f,7r},{f,|,0},{— f,f,0}  as  shown  in 
Figs.  l(b)-l(e).  The  ground-state  flux  pattern  depends  on 
the  relative  size  of  all  the  coupling  constants  in  Eq.  (1),  as 
seen  in  Fig.  2.  For  instance,  when  .p  <<C  { JA ,  Jv  ■  J  a  ■  -K'  1  • 
a  uniform  flux  state  of  {f ,f ,0}  or  —  |,0}  is 

favored. 

For  small  A,  the  ground  state  has  the  uniform  flux  config¬ 
uration  { j,j,0}  shown  in  Fig.  1(d)  and  all  its  excitations  are 
fully  gapped  for  a  system  with  periodic  boundary  conditions. 
The  total  flux  per  lattice  unit  cell  is  7r;  it  follows  that  the 
magnetic  unit  cell  must  be  doubled  by  a  gauge  choice. 
Note  that  even  though  a  specific  gauge  choice  it,j  with  this 
flux  configuration  would  “seemingly”  double  the  unit  cell 


FIG.  2.  (Color  online)  (a)  Ground-state  energy  per  site  of  (1)  as 
a  function  of  J5/J  with  JA  =  7/  =  Jy  =  =  J  for  the  uniform 

flux  configurations  shown  in  Figs.  1(b) — 1(e).  (b)  Ground-state  flux 
configuration  among  Figs.  1(b)- 1(e)  as  a  function  of  JA,Jv,Js  for 
JA  +  +  752  =  1  with  =  A  and  JA  =  JA.  The  coloring  scheme 

follows  that  of  the  legend  in  Fig.  2(a). 
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EXACT  CHIRAL  SPIN  LIQUID  WITH  STABLE  SPIN  . . . 


(a)  (b) 


FIG.  3.  (a)  Band  structure  on  a  cylindrical  geometry  for  JA  — 
J'A  =  1.0,  7v  =  7(,  =  0.8,  Js  =  0.  There  are  two  gapless  chiral 
Majorana  edge  states  (dotted  lines)  which  overlap  on  each  other, 
(b)  For  7A  =  1.0,  J'A  =  0.6,  7V  =  0.9,  7(,  =  0.5,  and  75  =  0.1,  the 
two  gapless  edges  states  separate.  These  ground  states  are  thus  CSLs 
with  a  spectrum  Chern  number  (±2)  and  the  vortices  obey  Abelian 
statistics. 

and  break  the  threefold  rotational  symmetry  of  the  lattice, 
the  physical  ground  state  obtained  after  the  projection  is 
translationally  and  threefold  rotationally  invariant  because  the 
gauge  choice  obeys  the  corresponding  projective  symmetry 
group  transformations  that  combines  a  physical  symmetry 
(here,  translation  or  threefold  rotation)  and  an  appropriate 
gauge  transformation.  ’4  Consequently,  no  symmetries  other 
than  the  time  reversal  are  broken  in  the  physical  ground  state. 

For  the  system  with  cylindrical  symmetry,  Figs.  3(a)  and 
3(b)  illustrate  the  existence  of  chiral  edge  modes,  indicating  a 
finite  spectral  Chern  number  of  ±2.  For  JA  =  J'A,  .7V  =  .A), 
and  J$  =  0,  the  c  and  cl  Majorana  fermions  decouple  and 
have  identical  spectrum.  It  follows  that  the  two  gapless 
chiral  Majorana  edge  states  overlap  as  shown  in  Fig.  3(a). 
Generically,  the  two  chiral  edge  modes  are  separated,  as  shown 
in  Fig.  3(b).  Moreover,  the  arguments  in  Refs.  15,  25  and  35 
for  spin  correlation  apply  here  as  well — the  spin  correlations 
vanish  beyond  nearest  neighbors,  establishing  a  CSL  phase 
with  a  finite  spectral  Chern  number  and  gapless  chiral  edge 
states  on  the  boundary  as  a  ground  state  of  (1).  Because  of 
the  even  spectral  Chern  number  (±2),  we  expect  the  vortex 
excitations  are  Abelian.  Nonetheless,  if  the  parameters  in  the 
Hamiltonian  doubles  the  unit  cell,  a  gapped  non- Abelian  chiral 
spin  liquid  can  also  be  obtained.36 

Gapless  spin  liquid  with  stable  spin  Fermi  surface.  For 
arbitrary  values  of  JA ,  Jv ,  7A ,  7(,  in  (1),  we  found  that  the 
system  is  always  gapped  as  long  as  we  kept  J$  —  0,  where 
Wp  =  —  1  for  all  p?6  By  increasing  ,/5,  we  found  that  the 
system  closes  its  gap  at  some  critical  values  of  .A  and  the 
ground-state  flux  configuration  changes  when  a  critical  point 
is  crossed  [see  Fig.  2(a)].  Note  that  the  energy  differences 
between  different  flux  states  are  surprising  small,  as  shown  in 
Fig.  2(a).  For  a  generic  7t  not  at  those  critical  points,  the  system 
stays  gapped,  which  seems  to  make  the  search  for  a  spin  Fermi 
surface  difficult  for  the  Hamiltonian  (1).  Fortunately,  while 
changing  7s,  a  particular  flux  configuration  can  be  always 
favored  by  adding  a  pure  flux  term  that  does  not  destroy  the 
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FIG.  4.  The  Fermi  surface  (solid  line)  for  the  flux  configuration 
{|,|,7r}and{7A. Tv, 75}  =  {1.0, 0.3, 0.8, 0.5, 1.4}.  Thedashed 
hexagon  is  the  Brillouin  zone  boundary.  Note  that  there  is  only  one 
Fermi  pocket  for  this  set  of  parameters  and  the  three  pockets  shown 
are  related  by  a  reciprocal  lattice  vector  and  are  thus  equivalent. 

exact  solvability  of  the  model.  Specifically,  to  stabilize  a  flux 
configuration,  we  add  the  following  term  to  Eq.  (1), 

H¥S  =  -aJ2w o-  (5) 

O  (A,  V) 

The  required  value  of  a  and  ft  to  stabilize  a  desired  flux 
pattern  depends  on  75,  but  is  typically  rather  small1617  as  the 
energy  differences  between  different  flux  states  are  small — see 
Fig.  2(a). 

For  appropriate  positive  a  and  ft,  the  ground  state  with 
the  flux  configuration  r]  (or  {—  f,7r})  is  favored. 

For  this  flux  configuration,  we  obtain  a  finite  Fermi  surface 
as  shown  in  Fig.  4(a)  for  7a  =  1.0,  7v  =  0.3,  JA  =  0.8, 
7y  =  0.5,  and  7,  =  1.4.  The  dashed  line  of  the  hexagon  is 
the  Brillouin  zone  boundary.  Note  that  there  is  only  one 
Fermi  pocket  around  one  inequivalent  corner  for  this  set  of 
parameters.  The  Fermi  surface  has  a  C 3  rotational  symmetry 
around  the  zone  corner  that  is  expected  from  the  C 3  symmetry 
of  the  model.  Varying  the  parameters  in  Eq.  (1)  by  a  small 
amount  will  only  change  the  shape  and  size  of  the  Fermi 
surface  but  the  Fermi  surface  itself  is  robust. 

An  important  question  to  ask  about  a  Fermi  surface  is 
whether  it  is  stable  against  any  weak  perturbations.  In  2D, 
a  generic  Fermi  surface  cannot  be  fully  gapped  by  a  weak 
“density  wave”  due  to  the  lack  of  perfect  nesting.  However,  a 
generic  Fermi  surface  can  still  be  fully  gapped  by  a  “pairing” 
term  with  momentum  Q  if  the  spectrum  at  k  and  —k  +  Q 
are  degenerate  for  any  k  and  some  fixed  Q.  In  general,  time 
reversal,  rr  rotation,  or  lattice  inversion  symmetry  enforces 
degeneracy  at  k  and  —  k  +  Q  and  allows  a  putative  Fermi 
surface  to  be  gapped  by  an  infinitesimal  “pairing”  term  with 
momentum  Q.  Note  that  a  nonzero  Q  is  encountered  in  the 
physical  state  but  no  gauge  choice  possesses  one  of  those 
symmetries.  If  the  ground  state  has  no  time  reversal,  n  rotation, 
or  inversion  symmetry,  a  putative  Fermi  surface  will  be  stable 
against  any  weak  perturbations.  In  the  current  model,  the 
7t  rotation  and  inversion  symmetries  are  explicitly  broken 
by  7a  7^  7y  or  J'A  f  7(,;  the  TRS  is  always  spontaneously 
broken.  It  follows  that  the  spin  Fermi  surface  found  in  our 
model  is  stable  against  any  weak  perturbations.  This  is  in 
contrast  with  the  spin  Fermi  surface  found  in  Refs.  15-17  that 
can  be  gapped  by  some  sort  of  weak  “nesting”  or  “pairing” 
terms.  Lastly,  this  state  has  recently  been  shown  to  exhibit  bond 
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energy  correlations  that  have  a  l/|r|3  power-law  behavior37 
despite  having  ultrashort  spin  correlations. 

Concluding  remarks.  We  studied  an  exactly  solvable  model 
of  spins  on  the  kagome  lattice  with  an  odd  number  of 
electrons  per  unit  cell  (making  it  a  “true”  Mott  insulator) 
that  realizes  a  gapped  chiral  spin  liquid  with  two  gapless 
chiral  Majorana  edge  modes  and  a  gapless  phase  with  a 
finite  spin  Fermi  surface  that  is  stable  against  any  weak 
perturbations.  While  in  terms  of  the  spin  variables  the  model 
may  seem  highly  anisotropic  and  therefore  somewhat  artificial. 
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Kitaev-type  models  have  been  derived  as  effective  low-energy 
theories  from  more  “realistic”  Hamiltonians  with  strong  spin- 
orbit  coupling38  and  can  be  readily  engineered  in  optical 
lattices.39,40 
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